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The model neutral exosphere with a uniformly rotating exobase is gener- 
alized by allowing variations in exobase density and temperature which charac- 
terize the thermosphere just below the base. The corresponding velocity 
distribution function, satisfying the collisionless Boltzmann equation, is con- 
structed and used to form a general expression for the velocity moments. 
Resulting density profiles of rotating exospheres with nonuniform densities and 
temperatures on the exobase are compared with corresponding nonrotating exo- 
spheres. Density enhancements due to rotation are found to be greatest above 
regions of exobase density or temperature minima. Equatorial density enhance- 
ments of terrestrial hydrogen, resulting from rotation, are estimated to be 15 
to 17 percent at altitudes of 10 to 20 earth radii. Corresponding increases in 
terrestrial helium are 30 to 50 percent on the equator at altitudes of 0.7 to 1 
earth radii even when there is a polar density bulge in the barosphere. 
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I. INTRO DUCTION 

The neutral planetary exosphere model of Opik and Singer 1 2 consists of a 
spherical base, or exobase, of radius R, through which all particles populating 
the collisionless exosphere emerge with a Maxwellian velocity distribution, 
characterized by uniform density and temperature and no rotation (no azimuthal 
bulk motion). On this model, Opik and Singer obtained an integral expression for 
the density profile which was subsequently confirmed by others , 3,4 through inde- 
pendent means, and placed in closed form. In a recent paper , 5 hereafter re- 
ferred to as I, Hagenbuch and iiartle extended the model by permitting the exo- 
base to rotate uniformly with angular velocity a. . It was found in I that important 
enhancements in exospheric density and escape flux on the equator (plane normal 
to rotation axis) can be affected by the inclusion of rotation. Burke 6 also noted 
increasing escape flux with rotation and discussed the possible astrophysical 
consequences. In addition to the rotational effects, important exospheric changes 
result when density and temperature variations are included on the exobase. The 
importance of such variations has long been recognized in the study of lateral 
flow in the terrestrial exosphere 7 ' 10 (flow from regions of high to regions of low 
exobase densities and/or temperatures). 

The purpose of this paper is to generalize the rotating exosphere model of I 
by allowing variations in exobase density, temperature, and bulk flow which 
characterize the thermosphere (barosphere) just below the exobase. We begin, 
in Sec. II, by constructing an explicit solution ol the collisionless Boltzmann 
equation from which we form a general expression for the velocity moments; 
e.g., density, flux, temperature, etc. The moments are in quadrature form and 
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must be treated numerically. In section III we present density profiles for 
selected conditions intended to bring out the principal effects resulting from a 
rotating exobase with density and temperature variations. Apart from that re- 
quired ior rotation, we ignore exobase bulk flow a? dom in all models to present. 
However, we discuss the possible importance oi such flow in the final section 
along with other results of the present model. 
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II. FORMJJt.i. ON OF MODEL 

The exobase surface is defined 2 by those points in the atmosphere where 
the lateral mean free path equals twice the density scale H; i.e., Iln = 1/2, in 
terms of the density n and gas-kinetic collision cross section . Then, due to 
the assumed variations in density and temperature, the exobase is a nonspherical 
surface with maximum and minimum radii. However, in this treatment we as- 
sume that the exobase is a sphere of radius H lying somewhere between these 
extremes. We also assume that the variations in base density and temperature 
are weak in the sense that the resulting difference between the maximum and 
minimum radii is small relative to R. 


We construct the velocity distribution function f in the inertial frame with a 
spherical coordinate system whose origin is at the center of the planet. In this 
case, the spatial coordinates are r, 6, and /, the radial, colatitude, and azi- 
muthal coordinates, respectively, and the corresponding Cartesian velocities 
are v r , v tU and v.,. At the exobase, r = R, the velocity distribution f (?, vf) of 
particles emerging into the exosphere is taken to correspond to local thermo- 
dynamic equilibrium to give the boundary condition 


f(R.v) 
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for v • R > 0 
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in terms of the particle mass m and Boltzmann’s constant k. The functional 
dependences of the particle density N 0\ + ), temperature T ( , /.), and bulk 
velocity U ( , ; ) characterizing the emergent distribution can be obtained from 
theoretical models or experimental data. Only simple forms for these functions 
are chosen in this work to illustrate the main features of this model. 


The exospheric distribution function t *\e seek satisfies the collisionless 
Boltzmann equation (in inertial frame) 
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in which the gravitational potential 

V = - MG/r, 

where M is the planetary mass and G the gravitational constant. The most 
general solution ol I.q. (2) consists of an arbitrary function of the five constants 
of particle* motion obtained by solving the corresponding characteristic equations 
to give 

L l - m (v ~ + Vg 2 + v 2 ) + mV, (3a) 
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The first three constants are the total energy, total angular momentum squared, 
and the z component of the angular momentum, respectively, while the remaining 
two constants relate the angular coordinates of a particle trajectory. 

In order to form the distribution function for the exosphere it is necessary 
to relate the point of observation (r, 0 % / ) of a particle of given velocity (v r , 

\ Q% v.) with the corresponding point on the exobase surface (R, 0, <t>) where the 
particle emerged; i.e., the source point. Such a relationship can be represented 
by the expressions 

rose = (1 - J J/J 2 ) cos [K + cos-' T,(R)] . (4a) 

cos *I> = cos [ L ♦ £ ( 1 - cos ? 0 ) ] , (4b) 

in term of the constants of motion (since 0 and <!> are functions of the constants 
of motion, they arc constant along a particle trajectory). For a point (v, r) in 
phase space we note that *.qs. (4) lead to four (E > 0) or eight (E < 0) points 
(©, <l>). A more specific representation, required in this work, is obtained by 
using Eqs. (3) and (4) and following the principles of elementary orbital me- 
chanics to give 

cos 0 i j cos 0 Y ( + e } P sin & Z x , (5a) 

oos 4\jk r { cos £ sin Y ( - [e j P cos^cos ; - f k ( 1 - P 2 ) sin />] Z.J/sinQ^, (5b) 
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where A (-1)' 1 for 1,2 and the four (eight) roots (ft, *h) • (W () ,‘!> ijk ) are 
Identified by the indices (i, j, k) through , sgn v r , sgn v , t and e. sgn v. . 
l he remaining functions are given by 

p I v $ ! (v/ ♦ v i)‘ , 

Y , r /(K) r 7 ( r ) f , t \(R) *( r), 

z , ^ ( r ) V ( R ) , i?(R) v(r), 

where 


X(r) [1 - ij a (r)]\ 


In addition to satisfying the boundary condition of Kq. (1), the distribution 
function w< seek only admits particles which have emerged from the exobase .and 
excludes those particles exceeding the escape speed with negative radial veloci- 
ties (i.e. t no source of particles at infinity), t he distribution function we pro- 
pose is 
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in which the unit step function 
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\ K is the potential at r R and the radial, latitudinal, and azimuthal components 

A 

of the velocity V are given by 

[i(E- jV2mR»-mV I ,)] M . 

J t mR sin 0 , 

respectively, and where the ± sign corresponds to sgn v n at the exobiuo. For 
simplicity, the subscripts (i, j, k) on ((-) , } , ( t’ ljk ) have been suppressed. The step 
function terms are analogous to those used by A a mod t and (Vse 4 and in I, in 
which case the step function multiplying the square bracket term admits only 
those trajectories having contact with the exobase while the term in square brackets 
removes all those particles exceeding the escape speed with negative radial veloc- 
ities. Altogether then, the distribution function of Eq. (6) admits only those trajec- 
tories allowed on the model and satisfies the boundary dondition of Eq. (1). 

In addition, it is clear that the kinetic equation (2) is satisfied since the distri- 
bution is a function of the constants of the motion except for the term involving 
S(v r l. It can be shown, by direct substitution into Eq. (2), that this term is an 
admissible solution (as shown similarly by Aamodt and Case). 

The velocity moments of the distribution are formed by multiplying f of 
Eq. (6) by the appropriate velocities and integrating over all velocity space. We 
simplify the analysis by normalizing the exobase density and temperature dis- 
tributions through 


N<a <i>) * n 0 v(a o>) , T<a <d> 


t 0 r(a«j>), 
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where n o ami l’ 0 are values of density and temperature at sHr -ted points on the 
exobasc . As mentioned above, we consider a corotating exobase of angular 
velocity "T directed along the polar axis and Ignore the effects of other possible 
hulk flows just below the exobuse by setting 

(l f r . U 9 . U A ) (0, 0, 1 ' K sin ft) . 

It is worth mentioning at thin ;>olni that difterentiai rotation couk be included 

here by letting < (^,‘1’). In addition, we make the spatial coordinate trans- 

formation 

r - K f 

and the velocity coordmat* transformation 


v 


r >" x * y“'*y. v r H z, 


in terms of the constant 


m 2 k T 0 . 


Then, the general velocity moment ot order (p, q, s), normalized to t^e h«i so 
density n 0 , is given by 
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correspond to contributions from the bounu and escaping trajectories, respec- 
tively. The limits of integration are given in terms of 

mMG / kT 0 K - 2> V R , 

a measure of the gravitational energy relative to the thermal energy at the base, 
and 

a = a/g, b = l/(f 2 - i), d a/(£ + 1), 

by 

P, = (P 2 2 - Z 2 )H, P 2 = (bx 2 + d)^, 

P 3 - [a(l - p 4 = (P 5 J - z>)\ 


P 5 = (a-x 1 )*, P 6 = a*. 


The integrand function 
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where 


,3 - (ma 2 RV2kT 0 )^ >*• . R 


is a measure of the rotational kinetic energy relative to the thermal energy at 
the base, and i/, lk *'(« ij ,< | > ijk ). and i, jk t(H, jV ). 


For completeness, it is worth mentioning at this point that the momenta of 
Eqs. (7) reduce to those obtaim d in I wh n the base U mperutui c* m taken to be 
uniform and the base density is taken to correspond to a rotating isothermal 
thermosphere (see Fq (8) below). In addition, for tht ease of a nonr Hating exo- 
sphere, 0, it can be shown through application of Tdouville's theorem and 
appropriate coordinate transformatlo: ; that Eqs. (?'■ reduce to the density* and 
radial flux 0 at the bas< , r R, as derived by McAffee. And finally, taking /? = 0 
and i ijk i, vve obtain the density distribution! of Opik and Singer. 1 ' 4 


The model considered here and in I is the usual minimum dens ty model ** 5 
in the sense that we ignore the possibility of bound orbiting particles (sate llite 
particles) which do not intersect the base (see rets. 2, i 1-13 for discussion 
of possible importance). The problem of determining the distribution of such 
particles is beyond the scope of this work (solution of kinetic equation which 
properly accounts for the production and loss of satellite p;irticles through rare 
collision events and photoionization in the exosphere). One may include the con- 
tribution of satellite particles in an approximate way by simply adding a Max- 
wellian distribution f^, truncated by [ 1 - S(E - J 2 /2mK ' - mV R )J S(-K), to 
Eq. (G). The temperature mid density characterizing f M could be taken as, say, 
those corresponding to the average values at the exobase. 
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IK. NUMER ICAL EXAMPLES 

A. Model I 

For purposes of brief review and comparison below, consider the model of 
I corresponding to a corotating, isothermal thermosphere where the exobase 
density distribution was given by 

N n Q exp (mco 2 R 2 sin 2 (-) t . /2kT 0 ) , (8) 

exhibiting a density bulge at the equator. We note that this density, given in 
terms of the present model, leads to an exospheric distribution functiDn f which 
is independent of ((-V ( , «l> )k ). Consequently, the analysis in I was greatly simpli- 
fied since detailed accounting of trajectories through such equations as (4) and 
(5) was unnecessary. In this case, the resulting profiles for density n and 
radial flux n v r •> simplify to a single quadrature over v 0 while n < > and 

n v, > vanish, (i.e., no lateral flow in exosphere). 

To illustrate some of the main features of the model in I, consider the 
dashed lines of Fig. 1 where the normalized density n/n 0 , at given colatitudes 9 , 
is plotted versus for the case i = 10 and 3 = 1 (these values have been selected 
primarily for clarity of illustration and are used throughout this section). As 
shown in I, the density distribution along the polar axis ( ;) = 0) is identical to the 
corresponding distribution of a nonrotating exosphere. We note, for a fixed 
radius, that the density increases from pole to equator. In other words, there is 
a density enhancement at all points in the exosphere for colatitudes 0 < 6 < n. 

As noted in I, this variation results since particles are ejected with increasing 
azimuthal velocities for increasing colatitude so that the average altitude a 
particle attains in the equatorial region is greater than in the polar region, 
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leading to a higher density at the equator than the pole. In what follows, we 
designate this phenomena as the "centrifugal effect.” 

B. Uniform Density and Temperature 

The conceptually simpler case of a corotating exobase with uniform density 
and temperature ( |jk ( ik 1) is computationally more difficult since f is a 
function of H ( and Eq. (7) is required for the moments. Such an example of the 
density profile for a 10 and 1 is shown by the solid lines in Fig. 1. The 
dashed line, designated by O' (the same as model I), is the corresponding 
density profile for the nonrotating case, valid for all latitudes. As in the model 
of I with rotation, the density increases from pole to equator for a given radius. 
Again, we note a density enhancement in the equatorial region which is primarily 
due to the centrifugal effect. On the other hand, the values of density along the 
polar axis are lower than corresponding values of the nonrotating case. In fact 
the density is lower for a range of colatitudes about the poles. 

This density reduction is consistent with the fact that there are now lateral 
winds in the exosphere flowing from the polar regions to the equatorial regions. 
Such lateral flow is expected in this case when a comparison is made with the 
isothermal model of I where no lateral flow occurs. In the latter model, the exo- 
base density, given by Eq. (8), increases from pole to equator in a manner which 
is consistent with no lateral flow. Then, in the former case, there exists an 
"effective" density depletion in the equatorial region of the base, relative to the 
isothermal case of I, leading to flow' from the polar region of effectively higher 
base density to the equatorial region of effectively lower base density. 
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C. Nonuniform Density 

Taking the simple uniform model of subsection B as a point of departure, 
consider a rotating exosphere characterized by an exobase with uniform tem- 
perature (t. . 1) and variable density. For simplicity, we choose an exobase 

density distribution which is rotationally symmetric about the x-axis ( / - 0, 

/2) with a sinusoidal density variation given by v ijk = N/n 0 = [ 1 - (1 - t N ) 

• (1 + < N )" 1 sin (•) i ( cos<l' jk ] , which results in a density minimum on the equator 
at 4> = 0 and a density maximum at 4> - v % The ratio < N N m , x /N m t n , in terms 
of the maximum and minimum exobase densities. 

Density profiles for a model exosphere of this kind, with - 1/3 (chosen 
for clarity), are shown in Figs. 2 and 3 for f> - v and = 0, respectively. In 
these figures n - 10 and the solid lines are for a rotating planet with /3 1 while, 

for comparison, the dashed lines are for the corresponding nonrotatir.g planet. 

In Fig. 2, note that the equatorial density is greater than the corresponding 
polar density at all altitudes above the exobase for both rotating and nonrotating 
models. This equatorial enhancement is in part a reflection of the exobase 
density maximum lying on the equator at this longitude. The ratio oi equatorial 

to polar density is more pronounced in the rotating exosphere due to a com- 
bination of the* centrifugal effect and lateral flow from regions of high to low 
exobase densities (actual and effective). 

For the nonrotating case of Fig. 3, we note that the polar density is greater 
than the corresponding equatorial density at all altitudes, reflecting the base 
density minimum on the equator at this longitude. In contrast, the opposite re- 
lationship obtains fer the rotating exosphere. That is, the equatorial density is 



1G 


greater than the* polar density at all altitudes above a crossover point x ~ 1.1. 

In addition to the centrifugal effect, this density increase on the equator is due 
to the "rotationally enhanced” lateral transport of particles from exobase regions 
of higher densities. In fact, such a "mixing” of particles tends to equalize the 
density profiles around the planet on a given latitude as one can note in com- 
paring the equatorial distributions of Figs. 2 and 3. 

To further elucidate this rotationally enhanced mixing, consider Fig. 4, 
where a sketch of flux field lines n v > are shown in the equatorial plane. For 
clarity, the relative number of field lines is not to scale since the flux decreases 
very rapidly with increasing . We note the strong connection of field lines from 
one hemisphere to the other where a portion of the field lines from the high den- 
sity base region feed back into regions of lower density while the remaining 
portion, corresponding to the escaping flux, form spirals which become radial 
as * 1 . The dashed lines shown indicate the approximate locus of points of 
maximum and minimum values of both flux and density for a given radial posi- 
tion. In this connection, we note that at about - 3, the density (;ind flux) maxi- 

mum is closer to the exobase density minimum than the maximum. This implies 
that the density at f - 3 should be greater above the base density minimum than 
the maximum; this is borne out by comparison of Figs. 2 and 3. 

D. Nonuniform Temperature 

We consider here the contrasting case of uniform base density (^ ljk = 
t N - 1) and a rotationally symmetric base temperature distribution r t jk 
T/T n = [l + (1 - t T ) (1 + e )~ 1 • sin 0 cos <t> . in terms of e = T /T , 
the ratio of the minimum to maximum base temperatures. Then, for the non- 
rotating model shown in Fig. 5 (6), with T 1/2, the equatorial density is 
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greater (less) than the corresponding polar density at all points above the base, 
this time reflecting the base temperature maximum (minimum) on the equator 
at this longitude. The same (opposite) relationship obtains in the rotating exo- 
sphere of Fig. 5 ((>). In these cases, the temperature maximum and minimum 
play similar roles to the density maximum and minimum of the previous ex- 
amples and the basic reasons for the differences between the rotating and non- 
rotating models are analogous. 


/ 
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IV. DI SCUSSION AND CON CH SIONS 

The exobase distributions in de nsity and temperature chosen for the numeri- 
cal examples of section III correspond to the common idealization in which a 
density minimum and temperature maximum occur on the equator /2 at 
the subsolar point , - 0 (approximate region of maximum solar extreme ultraviolet 
photoionization and heating) while a density maximum and temperature minimum 
occur at the anti-solar point (approximate region ot maximum recombination and 
cooling). Then, to maximize the effects oi rotation, the angular velocity * was 
directed along the polar axis. For purposes of illustration, the effects of vari- 
ations in exobase density and temperature were considered separately but, of 
course, occur together in nature. 

In the preceding examples we only considered the values i = 10 and ft - 1 , 
since the general effects of variation of these parameters, as examined in I, still 
hold for the present nonuniform models. For example, a primary result 
noted in I is that the ratio of the equatorial density to the polar density, at a 
given altitude and longitude, increases when , , and are increased separately 

or together. In addition, this ratio is greater than one for greater than the 
crossover point x . In this connection, we point out that it can be shown that 
- 1 as ft and/or a become large. 

On comparing Figs. 2, 3, 5, and 6, one can note that the region of greatest 
exospheric density increase from the nonrotating to the rotating model occurs 
above points of base density or temperature minima. Further comparison with 
the uniform model of Fig. 1 shows that the corresponding rotational density 
enhancement is less than that of the nonuniform models (above base density and 
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temperature minima; . Altogether then, the effects of rotation appear to be more 
important in nonuniform exospheres. 

With the model developed here wo can now estimate the equatorial density 
enhancement (region where most prominent) of terrestrial hydrogen and helium 
due to rotation. For this purpose we choose the distributions , , 1+0.28 

• cos 2>5 (Z. jk 2) and p. )k ; 1 — <l/3) cosZ ijk (here cos Z 1Jk «ainH i cos<t> ik ) 
used by Wallace, et al. M to lit the observed hydrogen density from Lyman i 
data. In addition, we assume that the exobase is corotating with , 1.02 x 10 4 

sec' 1 , the observed atmospheric (super) rotation rate 15 near the exobase for the 
range 30° £ (■ i 150 ; i.e., 1.4 times the rotation rate of Earth. In this case, 
with R = 0880 km and T 0 805°K, we obtain a 15 to 17 percent increase in hy- 
drogen density above that of the nonrotating model for the range of 10 to 20 earth 
radii. 

The corresponding density enhancement in terrestrial helium is 40 to 80 
percent over the range 1.7 to 2 planetc centric radii, the region where helium 

9 

becomes a minor constituent. However, recent observations indicate that there 
is a helium density bulge over the winter pole. 16 In this case, distributions of 
the form v t) = 1 + .82 cos (4 , and r. * 1 may be more realistic. On this 
basis, taking T 0 1000 K, w'e still obtain a rather large density enhancement 
of 30 to 50 percent on the equator for c of 1.7 to 2. 

The numerical examples considered in this paper assume that appropriate 
sources and sinks exist in the thermosphere to allow the resulting flow of par- 
ticles across the exobase. The flow is limited by the maximum rate at which 
the species in question can diffuse in the thermosphere. No quantitative 
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limitations can be determined for the examples of Sec. ill until a planet and 
species are specified. However, for the terrestrial hydrogen ;ind helium ex- 
amples of this section, the maximum allowable iliffusivc fluxes, as estimated by 
McAfee, are about 8.9 10 7 cm ’sec' 1 and 5.1 • 10* cm’*’ see 1 , respectively. 

In these cases, we obtain maximum allowable exobase densities of 5 x 10 4 cm 3 
(somewhat lower than the expected value 10' cm 3 ) ano > 10* cm 3 (expected 

value 1() 6 cm 3 ) for hydrogen and helium. 

Exosphere models of the kind presented here assume a Maxwellian velocity 
distribution at the exobase. This assumption is self-consistent with respect to 
the distribution of bound particles, E 0, when the exobase is uniform. On the 
other hand, as is well known, the distribution has a "hole” corresponding to the 
absence of those particles with E 0 and v t 0, leading to a flux (and density) 
discontinuity across the exobase, the Jeans escape flux. With a nonuniform exo- 
base, even the distribution of bound returning particles, I 0 and v r 0, is not 
Maxwellian in form, leading to further discontinuity in density and flux across 
Ihc exobase. But, when variations in density anti temperature occur just below 
the exob.isi , flow is expected in the thermosphere and should be reflected in the 
exobase distribution. The possibility for such flow was included in the general 
distribution of Kq. (G). 

In the context of this model, one cannot expect to remove the density and 
flux discontinuity corresponding to those particles in the range E ■ 0 due to the 
complete absence of such particles with v r 0. However, since bound returning 
particles (E 0, v r 0) are present, it is logical to expect that it should be 
possible to remove or reduce the discontinuities in flux and density of the class 


of particles with F 0 by balancing the tlux just below the base {bound particles 
only) with the flux n 0 (‘A, 00 r + rt ((10 ■ + rt 00 , / ) just above (conservation of 
flux). For example, when the density and temperature are given at the base, 
the balancing can be accomplished through successive approximations to 

U (H, <l>). To test this procedure and get an indication of the effects of including 

• • 

bulk velocity U, we performed approximate hand calculations to balance the flux 
across the base (10 percent accuracy) for the examples of Fig. 2. We found that 
the density discontinuity was reduced and the density profiles were essentially 
unchanged for y >10" 13 (over the range considered 1 - _ 10). When y was 

reduced to 10" 14 , we found both increases ;ind decreases in the density by about 
10 percent at 10 ; such changes being smaller for f 10 and vanish as ; - 1. 
In addition, the maximum outward radial bulk velocities U r were 5 x 10 4 cm/sec 
and 1 x 10 s cm/sec when / was lO" 13 and 10~ 14 , respectively. Also note that, 
when U r > 0 on this model, the modified Jeans escape flux n 0 J l00 is greater 
than the corresponding value when U f 0, which may result in important in • 
creases in estimates of planetary loss rates. 

The degree to which the base can be nonuniform on this model and still lead 
to meaningful solutions is unknown. As done previously for the uniform 
model, 17 * 19 this problem can be resolved by a full kinetic treatment of the exo- 
base region, taking proper account of the transition from the collision dominated 
thermosphere to the collisionless exosphere. It is clear, of course, that the 
accuracy of the model improves in the limit of a uniform exobase, in which case 
U 0 (distribution of bound particles becomes Maxwellian at base). 
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FIGURE CAPTIONS 


Fig. 1. 1 he normalized density n/n 0 at colatitudes 5 0°, 45 , 90° versus the 

normalized radius ft - r/R for values n - mMG/kT 0 R = 10 and ft - 
(m . 2 R 2 /2kT 0 )‘* = 1.0. The dashed lines are the isothermal model of I and 
the solid lines are for the model with uniform base density and temperature. 

Fig. 2. The normalized density n/n 0 at colatitudes Q = 0°, 90°, longitude b = 180° 
versus the normalized radius f = r/R for values \ - mMG/kT 0 R 10, 

(m - 2 R 2 /2kT 0 )'* = 0 (dashed lines), 1.0 (solid lines), and < N = N n t / 

N = 1/3. 

Fig. 3. The normalized density n/n () at colatitudes & - O’, 90°, longitude b = 0° 
versus the normalized radius ft = r/R for values i = mMG/kT 0 R = 10, 
(mo 2 R 2 /2kT 0 )^ = 0 (dashed lines), 1.0 (solid lines), and < N = N min / 

K.. • 1/3. 

Fig. 4. Sketch of field lines n < v > in the equatorial plane for counterclockwise 
rotation. The maximum exobase density is indicated by II and the minimum 
by L. The dashed lines to the right (left) are the approximate locus of both 
density and flux maxima (minima) for a given altitude. 

Fig. 5. The normalized density n/n 0 at colatitudes 0 = 0", 90°, longitude b - 0° 
versus the normalized radius ft = r/R for values a = mMG/kT Q R = 10, 

? = (ma) 2 R 2 /2kT 0 ) 4 = 0 (dashed lines), 1.0 (solid lines), and T = T min / 

T,... = 1/2. 


. f). The normalized density n/n (1 at colatitudes = 0% 90% longitude / 180° 

versus the normalized radius 1 r/K for values » - mMG/kT () R 10, 

(m- 2 R 2 /2kT 0 ) • 0 (dashed lines), 1.0 (solid lines), and < T T n . in / 











